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Four dimensional simplicial gravity has been 
studied by means of Monte Carlo simulations for 
some time Q, the main outcome of the studies 
being that the model undergoes a discontinuous 
phase transition between an elongated and a 
crumpled phase when one changes the curvature 
(Newton) coupling. In the crumpled phase there 
are singular vertices growing extensively with the 
volume of the system |3| giving an infinite Haus- 
dorff dimension, whereas the elongated phase has 
a Hausdorff dimension equal to two. This phase 
has all properties of a branched-polymer phase 
H). We have postulated]^ that this behaviour 
is a manifestation of the constrained-mean-field 
scenario as realised in the Branched Polymer]^ 
(BP) or Balls-in-Boxes model[0. The models of 
jsj-Q share all the features of 4D simplicial grav- 
ity except that they exhibit a continuous phase 
transition. We note here that this defect can be 
remedied by a suitable choice of ensemble. 

The partition function of the Balls-in-Boxes 
model [fzl is 



Z{M, N) 



piqi)---piqM)6q,+...+q,,^N, (1) 



which describes weighted partitions of N balls in 
M boxes. This model can be solved[0 in the limit 
of an infinite number of boxes and fixed density 
of balls per box : M — > cxd and p = N/M = 
const. In this limit the partition function can be 
expressed in terms of the free energy density per 
box /(p) 



Z{M, N) 



(2) 



Kronecker delta function one finds by the steepest 
descent method that 

f{p)^fi4p)p + K{fi4p)) (3) 
where if is a generating function given by 



(4) 



if(Ai) = log^p(<z)e-^« 

9=1 

and p.*{p) is a solution of the saddle point equa- 
tion 



p + K'{p^) = 



(5) 



For a suitable choice of the weights p{q) the sys- 
tem displays a two phase structure with a critical 
density pcr- 

When p approaches pcr from below /x* ap- 
proaches Per from above. When p is larger than 
Per , P- sticks to the critical value pcr and the free 
energy is linear in p 



f{p) ^ PcrP + Kc 



(6) 



By introducing the integral representation of the 



where Hcr = K{pcr)- The change of regimes 
P < Per (^ to p > Per (||) corresponds to the 
phase transition. To understand the physical na- 
ture of the transition it is convenient to consider 
the dressed one-box probability, defined as the 
probability that a particular box contains q balls. 
In the large M limit the saddle point equation 
gives 

r e-^(''-)p(q)e-'^*'!' for p < p^r 

<q) = <^ (7) 
[ e '^"'p{q)e for p > pcr 

The approach of the dressed probability to the 
limiting form (^ is not uniform. In particular 
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the average, (q) — J2q QT^il): does not give p for 
p > Per as it should. One can easily correct for 
this by adding a "surplus anomaly" term to n for 
finite M in the phase above pcr [@| 



M{p - Per))- 



(8) 



In effect, M — 1 boxes keep the critical form of 
the distribution and one box takes over the sur- 
plus of balls. One can easily check by performing 
directly finite size computations that this situa- 
tion is indeed realized in the model. 

The appearance of the surplus anomaly is 
a condensation phenomenon similar to Bose- 
Einstein condensation and the transition in the 
Kac-Berlin spherical model We call the phase 
P < Per fluid and p > pcr condensed. The sys- 
tem may enter the condensed phase either by 
changing density or by modifying the weights p. 
For instance, the one parameter family of weights 



= 9 9 > 1 gives the critical line 



C(/3 - 1) 
C(/3) 



(9) 



so one can change the phase either by varying (3 
for fixed p (as in case of branched polymers where 
/9 = 2), or by fixing (3 and varying the density p. 
The transitions in p and /3 in this variant of the 
model are continuous, unlike the simulations of 
simplicial gravity. 

It is also possible to consider ensembles with 
varying density. The simplest candidate is an en- 
semble with a chemical potential coupled to the 
total number of halls^ which can be treated as a 
box in contact with a reservoir of balls. However, 
this leads to a totally decoupled system which cor- 
responds to M copies of the urn- model For 
P > Per the average number of balls in the urn di- 
verges. More interesting in the context of simpli- 
cial gravity is the model where we keep the num- 
ber of balls fixed and vary the number of boxes 
[^,|ri[. This gives direct analog of the ensemble 
used in the simplicial gravity simulations 



Z{n,N) ^^Z{M,N)e' 



M 



(10) 



It is now more natural to consider curvature r — 
1/p instead of density. The partition function 
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Figure 1. The average curvature as a function of 
K. The heavy line corresponds to thermodynam- 
ical limit. {(3 = 5/2) 



( p^ ) can be rewritten as 
1 

Z{k,N)kN j dre^U(.r)+^r) ^^^^ 

l/N 

and the saddle point equation for this integral is 

«. + /'(n)-0. (12) 

For K > Her this equation reduces to ft = 
K{psp{r*)) which has a unique solution for r^,. 
The value of r, (the centre of a gaussian dis- 
tribution) is the average curvature in the limit 
N —>■ oo. This situation continues as long as 

For K < Her the saddle point equation ( p^ 
has no solution and therefore the integrand is 
no longer gaussian but a monotonic function of 
r. In particular, for r < r^^ it is exponential, 
exp N{k — Kcr)r, and for large N only this expo- 
nential part contributes in the integral (|l^) giving 
(r) ~ l/N. 

The average curvature (r) is shown as a func- 
tion of K in Figure |^. The bold line is a limiting 
curve for N = oo. For k > Ker it is the solu- 
tion of the saddle point equation (p^). It stops 
at rer and falls to zero. In the neighbourhood of 
the critical point the curves are steepest. This 
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Figure 2. Distribution of r near the phase tran- 
sition for systems with 512 (k = —0.32184) and 
1024 balls (k ^ -0.31910) {(3 ==5/2). 



part of the curves corresponds to the pseudocrit- 
ical region where the two phases coexist. One 
expects a double peak histogram for r : one peak 
near the maximum of the gaussian phase and the 
other near the kinematic limit 1/A^. In Figure |^ 
we show the distributions of r for two different 
sizes N, showing the coexistence of two phases 
that is characteristic of first order transitions. 

To summarize, the Balls-in-Boxes model de- 
scribes well basic features of simplicial gravity 
simulations such as the appearance of the sin- 
gular vertices and the mother universe With 
the appropriate choice of ensemble one obtains 
a first order transition too. We note that argu- 
ments 14,1^] have recently been given that the 
bare weights p{q) in simplicial gravity can be ap- 
proximated by p{q) ^ e^'^e^"''^ , which gives a 
similar behaviour to the power law weights dis- 
cussed here. 

This work was partially supported by the KBN 
grant 2P03B04412. 
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